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Abstract 

The differential equation method is applied to evaluate analytically two-loop 
vertex Feynman diagrams. Three on-shell infrared divergent planar two-loop 
diagrams with zero thresholds contributing to the processes Z — > bb (for zero b 
mass) and/or H —>■ gg are calculated in order to demonstrate a new application 
of this method. 



1 Introduction. 



In the two-loop order of the Standard Model large effects are due to the exchange of heavy 
particles, like H and t. The present investigation is motivated by the need to perform 
accordingly high precision calculations of the decay rates Z — > bb and H — *> gg. We 
demonstrate a new application of the differential equation method (DEM) to calculate 
analytically several scalar vertex diagrams for the above decays. Intense study of the 
relevant diagrams is in progress (see 0, H) : semianalytically, using the technique 
of asymptotic expansion of Feynman diagrams 0, |^] and subsequent use of conformal 
mapping and summation by Pade approximation. We would like to note, however, that 
the technique of [|l], works best for diagrams with high thresholds. In the case of several 
massless virtual particles and in particular if the diagrams under consideration have zero 
thresholds, additional contributions (apart from the "naive" diagram, see e.g. 0) are 
needed due to the underlying method of the large mass expansion. Moreover, in the above 
approach the convergence of the series w.r.t. the inverse masses is getting worse for large 
external momenta squared the fewer heavy masses are involved. Thus, to obtain analytic 
results is quite an important task in particular for those cases with many massless virtual 
particles. If only one non-zero mass is involved, in fact the DEM provides in a systematic 
manner these analytic results in a managable amount of work. These analytic results in 
general are given in terms of a one- dimensional integral representation from which it is 
possible to obtain the Taylor series expansion in the external momentum squared. The 
latter is considered as the most efficient approach for the final numerical evaluation due 
to references |T], |3| . 

In the next section we shall recall the basic formulae and the sence of the DEM f|-[|| 
(and |9j for short review), which is used for the calculations. In Sections 3 and 4 we shall 
demonstrate the basic steps in the calculation of the massless two-loop diagram (Case 10, 
according to the notation of [[| and Fig. 1) and one diagram with one non-zero mass 
(Case 7). Section 5 contains only the results for another diagram (Case 8). 



2 The Differential equation method 

The DEM is a method which allows one to obtain results for massive diagrams, evaluating 
only diagrams having an essentially simpler structure than the initial ones. The method is 
based on the rule of integration by part |I| . Throughout the article we use the following 



notation. Assuming dimensional regularization, all the calculations are performed in 
Euclidean momentum space of dimension D = 4 — 2s. The dotted and solid lines of any 
diagram correspond to massless and massive euclidean propagators [] 

(7f = ' (q 2 + m 2 ) a = 

a and m are called the index and mass of this line. Further (except when mentioned 
otherwise) all solid lines have the same mass m. Lines with index 1 and mass m are not 

lr The transition to Minkowsky space for a diagram / — > / is given by changing qo — > iqo, qj — > qj (j = 
1,2,3) (q and q are the external momenta of the diagram in Euclidean and Minkowsky space, respectively) . 
Hereafter all variables with " ~" denote Minkowsky space. For the diagrams under consideration, we will 
make these transitions at the final step of the calculation. 
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marked. Let us now consider the rules for our calculation. 



Rule 1. Massive tadpoles T(«i, a 2 ) are integrated due to the graphic identity (here- 
after we use the Euclidean integral measure d D k/ix D / 2 ) 



R(a 1 , a 2 ) 



where 



2^ ai +a 2 -D/2 

T(D/2 - oti)r(a a + a 2 - D/2) 



Y{a 2 )V{D/2) 
and T is the Euler T-function. 



Rule 2. For a triangle the following recurrence relation |10 
line with index «j has mass 777.3) 



is valid (here the 




(D — 2azi — a 2 — a 3 ) 



Oil 



+a 2 \ 



«2 + l /\ a 3 




"2 + l/\ Q 3 



Ol — 1 




-2m\ai 



—{ml + m 2 ] 




ai+l 

02 + 1 /\ a 3 



-1 Ql 




We stress the fact that the basic line of the triangle always plays a special role (in the 
following we call it "distinctive"). 

This Rule, i.e. the rule of integration by part, is obtained by multiplying the integrand 
of the original diagram by D = dk^/dk 11 (k is the momentum along the distinctive line) 
and using / dPk div() = for the regularized Feynman integrals. 

Rule 3. Massless loops L(q 2 ;ai,a 2 ) with the external momentum q are integrated 
due to the graphic identity 



"2 



where 



A(ai, a 2 ) 



a(ai)a(a 2 ) 



T(D/2 - a) 



a{a.\ + a 2 — D/2) ' r(a) 

Stricly speaking, Rule 3 is not independent of Rules 1 and 2 because it may be ob- 
tained (at least for even values of a± or a 2 ) in the limit m 2 — > from equation (1) below. 
Hovewer, it is very convenient to consider it as third additional rule. 

Consider a few specific examples. 



1. A simple loop with one massive line 

= h(q 2 ;a) 
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Applying Rule 2 {a,\ — 1, a 2 — a, «3 = 0) with the massless line as distinctive one, we 



get 



(D -2- a)h{q 2 ]a) 



a 



T(0, a + 1) - {q 2 + m 2 )/i(g 2 ; a + 1) 



The last diagram — ali(q 2 ; a + 1) on the r.h.s. is the derivative with respect to m 2 of the 
initial diagram Ii(q 2 ; a). Hence, the r.h.s. has the form 

aR(0, 1 + a) +!_ g/2 + (g 2 + m 2 )-^h(q 2 ; a), 

Thus we get a differential equation with respect to m? for the initial diagram. Its solution 
for Ii is 



g 2(a + l-D/2) /i(g 2. 



v 1 /■* ai?(0,l + a) , g 2 



Notice that sometimes a more convenient representation of the diagram I\ is given in the 
form 



aR(0,a + 1) 



ds 



l-D/2 



(1 _ s y-D/2 s a+l-D/2 q m 2 /s 



2. A simple loop with two massive lines 



/ 2 (g 2 ) 



In analogy with the previous consideration we get 



->3-L> 



R(Q,2) 



ds 



1 _ s )l/2 s 2-0/2 g 



2-D/2 



4m 2 /s 



It is clearly seen that Eqs.(l) and (2) allow us to reduce an Z-loop diagram to an 
(I — l)-loop diagram with one propagator having "mass" m 2 /s and Am 2 /s, respectively. 
Applications of these relations may be found in || 0. 

Notice that if some lines have index 1, Rule 2 leads to a differential equation for the 
initial diagram with simpler diagrams in the inhomogeneous term. If, however, all indices 
are different from 1, Rule 2 does not lead to a simplification of the inhomogeneous term. In 
this case a more complicated technique, for example the method of Feynman parameters, 
is needed. Using the method of Feynman parameters the relation between /-loop and 
(/ — l)-loop diagrams in general may be obtained in the form (as before corresponds 
to the mass ra,-): 



cti , mi 




T(a 1 + a 2 - D/2) f 1 ds 

(1 _ s \ai+l-D/2 s a 2 +l-D/2 q 



ai+a 2 -D/2 

1 2 
m l m 2 



(3) 



For «i = 1 and a 2 = a. and m\ = (3) reduces to (1). 
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3. Vertex-type diagrams and iV-point functions were studied in |7| and ||, respectively. 
Further simple examples are worked out in these Refs. Since the basic subject of our 
investigation is to calculate vertex diagrams in Z — > bb and H — > gg decays, we consider 
here as another example the one-loop vertex diagram with the kinematics q\ = q\ = 
0, (qi - g 2 ) 2 = q 2 - 

qi - fta 

/ \ = I 3 (q 2 ,m 2 ) 

q-i/ \qi 

Applying Rule 2 with the massive and one of the massless as distictive lines, we obtain, 
respectively 

A A 

(D - 4)/ 3 = 2L(q 2 , 1, 2) - 2m 2 / % \ -2m 2 / \ (4) 

A k 

A A 

(D - 4)/ 3 = T(l, 2) + T(2, 1) - q 2 / '\ -m 2 / \ (5) 

' v 

a) For the massless case, Eq.(4) leads to the result for Is(q 2 , 0) in the form 
2 1 ah c\ 1 1 T ( 1 + £ ) r2 ( X ~ £ ) 



W,0) = --A(l,2) 



e K ' V (1+e) e 2 ? 2(i+ e ) r(l-2e) 



b) For the massive case combining (4) and (5), we get the differential equation for 
Iz(q 2 , m 2 ) in the form 



2m 2 



(£>-4)J 3 = 2/ 3 + 2m 2 1 



m 2 \ d 



q 2 J dm? 



where 



h = L(g 2 ,l,2)-^(r(l,2) + T(2,l)) 



is a combination of the massless loop and massive tadpoles. 

Introducing the variable x = q 2 /m 2 and using the boundary condition: J 3 = for 

m 2 — > 00, we get 

h{q\m 2 ) 



2; 



dy y 



-2c 



lT(l + e) x 

e g 2 ( 1+£ ) (1 - x) £ Jo (l-y) 1 - 

i[c(2)-Li 2 (l-rr; 
q L 



r 2 (i-^) 

r(l - 2e) 



(6) 



where Li n is the polylogarithm and ((n) the Riemann ( - function. 
In Minkowsky space @ reads 



J 3 (g 2 ,m 2> 



C(2)-Li 2 (l + £) 
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where x = q 2 /m 2 . The appearence of an imaginary part for 2/3 from Li2(l + x) for x > 
is seen from the analytic properties of Li2 . 

The key idea of the DEM: using the rule of integration by part (Rule 2) with different 
distinctive lines, the results for massive Feynman diagrams are obtained as a differential 
equation with an inhomogeneous termQ containing simpler diagrams. Using this procedure 
repeatedly one can reduce the 2-loop diagrams to 1-loops. These are integrated either 
again with the help of Rule 2 (by analogy with example 3), decreasing step by step the 
number of lines, or using the method of Feynman parameters for more complicated cases. 

The main difference from the massless case is the necessity to integrate the final 
result with respect to the mass several times. The basic advantage compared to the 
usual methods (for example, Feynman parameter method) to calculate massive Feynman 
diagrams is the appearence of complicated functions in the final step only. Hence, these 
functions do not interfere the process of calculation. 



3 Two-loop vertex diagrams 



At two- loop level we consider here three diagrams (see Fig. 1) out of ten characterized 
in They have the kinematics: qf = q\ = 0, (q± — (72) 2 = q 2 . The evaluation of the 
massless diagram J\q {Case 10) and one diagram with a non-zero mass J7 {Case 7) will 
be fully presented. For another diagram (J 8 , i.e. Case 8) we qive only the final result. 
Applying Rule 2 with line 6 (see Fig. 1) as distinctive one, we have for the J-j diagram: 



{D-A)J 7 = 2 






2m z 




(7) 



Repeating the application of Rule 2 to the second diagram of Eq.(7) with the left 
vertical line as distinctive one, we get 



p-4) 




Case 10. For the massless case, Eqs.(7) and (8) lead to the result for J 10 in the 
following form 



-2e)J 10 = 2,4(1,2) 



2(1+6) 



• \ \/ -1/4(2,2) \J +t 

y / y 



(9) 



The first diagram on the r.h.s. of Eq.(9) was already calculated in the previous section, 
the third one may be done in an analog manner. The second diagram may be evaluated 

2 Sometimes (see Eq.(8) below, for example) the differential equation can be brougt to a simple arith- 
metic relation between the initial diagram and simpler diagrams in the inhomogeneous term. 
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by the method of Feynman parameters. The result for Jiq is known (see 
form: 



J 
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+ 



1 r(l + e) T\l-e) 
£ a ( g2 )2(i+ £ ) r(x _ 2e ) 

3r 2 (l-2e)r(l + 2e) 
4 r(l - 3e) 



r 2 (i- £ )r(i + g ) 

r(l-2e) 



and has the 
3r(l-2e)r(l + 2e)r(l-e) 



1 r 2 (l+e) 

4e 4 g 2(l+ £ ) 



1 + 4C(2)e 2 + 20C(3)e 3 + 22((4) £ 4 



r(l-3e)r(l + e) 



+ 0(e) 



(10) 



We consider our derivation particularly simple. In Minkowsky space the diag ram J \q is 
represented in the form (10) with the replacement q 2 — > —q 2 and the additional factor —1. 

Case 7. For the case with one non-zero mass, Eq.(7) leads to a differential equation 
for J7 since the sum of the last two diagrams on the r.h.s. is the derivative of the initial 
diagram J7 with respect to m 2 . Solving it, we have 



x dy 



y 



l+e 



ft 



'.ID 



where fj is the inhomogeneous term containing the two first diagrams of the r.h.s. of 
Eq.(7). These diagrams are effectively only one- loop vertex diagrams. Indeed, using for 
the second one Eq.(3) with m\ = to the first two diagrams of the r.h.s. of Eq.(8) and 
with m\ = m 2 , to the last one, we have 



l+e 



r(i + e) 

2e 



ds 



(s(l -a)) 



l+e 




-(l + e 




(12) 



Evaluating the one- loop diagrams on the r.h.s. of Eq.(12) by the method of Feynman 
parameters and calculating the one-loop parts of the first diagram on the r.h.s. of Eq.(7), 
we obtain the final result for the diagram J7 from Eq.(|ll]) (neglecting terms ~ 0(e)) 



J7 



r 2 (i + £ ) 

\2e 



r l rja \ 1 I/I \ 

/ Sr^ ? 1 "(l^) + -^ln 2 (l + ,)-ln(x).n(l + ,)) 



2(q 2 ) 2 (m 2 ) 2£ Jo 8(l-s) 
+ <S> { J\z) + - ln(x) ln 2 (l + z) + \ ln 2 (x) ln(l + z) 



(13) 



where z = xs(l — s), 



^J\z) = \n[s(l-s)]\n 2 (l + z) - Mn 3 (l + z) - l2S l>2 (-z) - 5 ln(l + z)U 2 (-z) 
- C(2)ln(l + «) 
and hereafter lO 
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Snb.(x) 



-1) 



n+k— 1 „i 



dy 

(n — l)\k\ Jo y 



\n n -'(y)\n k (l-xy), S n>1 (x) 



Li 



n+l 



X) 
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1. In Refs. Jl|, H, [| some diagrams contributing to the processes Z — > 66 and H ^ gg 
were evaluated semianalytically, i.e. the first coefficients of the Taylor series in q 2 /m 2 were 
calculated analytically (in terms of rational numbers and C(2), C(3)) and the diagram itself 
is reconstructed by means of conformal mapping and the Pade approximation technique. 
It turned out that this approach is particularly useful for the final numerical evaluation 
since it allows in general to calculate the diagram to higher precision than the numerical 
integration. Therefore we give below the Taylor series of (13). The importance of the 
obtained formula below is that arbitrary many Taylor coefficients can be obtained very 
easily while in Ref.|| their calculation needed much computertime. This means that with 
such a representation the diagram can be calculated numerically to even higher precision 
with almost no effort. 



T 2 (l + e) 



E 



x) n T 2 (n) 



(g 2 ) 2 (m 2 f ti r(2n + l 



ln(x) + Si(n — 1 



:So(n-l) 



—S 2 (n ~ X ) + AS ^ n ~ i) 5 ^ 2 ™) - C(2) - Hx)S 1 {n - 1) + X - ln 2 (x) 



where 



Si(n) 



n ^ 



2. In Minkowsky space the diagram Jj is represented in the form (13) with the 
replacement x — > —x and the additional factor —1. The imaginary part generates from 
the terms ln(— x) = ln(x) — in and ln 2 (— x) below the threshold at q 2 = Am 2 and from 
the functions S^^rr^l — s)) above. 

At threshold q 2 = Am 2 we can get a closed form for diagram J7. Splitting J 7 into real 
and imaginary part: Jj = Re J7 + ilm J7, we have 



Re J 7 



1 



7T 



Im J7 



T 2 (l + e) 



16(m 



3 <(2) - -(7C(3) + 6C(2)ln(2) 
511 



+ 6C(2) ln 2 (2) - 14C(3) ln(2) - —((4) + 16^,1 



r 2 (l + e) 
16(m 2 ) 2+2e 



<(2)-6C(2)ln(2) + 7C(3) 



where the function 

U31 = 



E • 

n>m>0 



n 3 m 



ic(2)ln 2 (2)-^ln 4 (2) + ic(4)-2Li 4 (i) 



was defined by Broadhurst in [13 
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Case 8. The scheme for calculating diagram J 8 is similar to the one of the previous 
subsection. Here we will present only the final result for this diagram. 



l.The result of the diagram J 8 may be given in the following form 

T 2 (l+e) 



where <5>f\x) = $^(ar) + and <S>£>(x) = $^(x) + &o>(x) 

and 



J 8 
>( 8 ), 



(q 2 ) 2 (m 2 ) 2e 



(14) 



■,(8), 



7(8) / 



i(8), 



$f(x) 
l>! 8) (x) 



C(2)-Li 2 (l-x), 

3Li 3 (l -x)- 35 , i, 2 (l - x) - 2C(2) ln(l - x), 



o 1 - y 



Li2(-y) 



+ C(2) 

= 4 f-^ 
Jo 1- 

- $S 8) (x)ln(l -x) 
or the expansion w.r.t. q 2 /m 2 : 



ln(l - x)Li 2 (-x) + -S 1;2 (x 2 ) - S 1>2 (x) - S 1>2 (-x) 
2C(3) - 751,3(1 - x) - 9Li 4 (l - x) + 11S 2 , 2 (1 - x) - ^Li 2 (l - x) 



ln 2 (l -x) + 6Li 2 (l -x) 



y 



-Li 3 (-y) + 4S lt2 {-y) + ln(l - y)Li 2 (-y) 



r 2 (i + £ ) 

2£ 



oo n 

E^ 



(^(m 2 )" ^ n 



i n 

72 



+ o 8) (n) - ln(x)0i 8) (n) + t l n 2 (x) [5^ + 



n 
1 



\\ 








7 


| 7 


n 


~ 6 



+ -ln 2 (x 

2 v ' 



In^x) 



where 



and 



# 8 >(n) = 35 2 -4K 2 + ^ + ^-C(2), 

n rr 

7 „ 9 55*i 7 _ 



q o2 _|_ ij c c >7 

0i 8) (n) = 9S 1 S 2 + 2S 3 -2X 3 -12X^-4^X2 + ^? - + 54 + " 



In 



n 2 n 3 



C(2) 75i + - - 2C(3) 



n-1 



n-l (i\k+l 

E 4- and tf 2>1 
1 ft 1 



_i\fc+i 
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2. In Minkowsky space the diagram J§ is represented in the form ([14]) with the 
replacement x —* —x and the additional factor —1. The imaginary part generates from 
the functions S n) k(l + x), when x > 0. 

The calculation of the diagram J 8 = Re J§ + ilm J 8 at threshold q 2 = m 2 leads to the 
following result 



Re Js 



— Im J 8 

7T 



+ 



+ 



r 2 (i + £ ) 

(m 2 ) 2+2e 
-(- 

r 2 (i + £ ) 

fm 2 ) 2+2e 



^C(2) 



C(3) - 15C(2) ln(2)) + ^C(2) ln 2 (2) - ^<(4) + y£/ 3 ,i 
1 



3M2) 



2£ 



53 



C(2) - ln'(2) + - ln<(2) - 14C(2) ln(2) + -C(3) 



4 Conclusion 

The differential equation method has been applied to demonstrate a new approach to 
calculate two-loop vertex-type diagrams with one non-zero mass. Analytic results have 
been obtained for several scalar two-loop diagrams of the processes Z bb and H — > gg, 
which will be relevant for tests of the Standard Model. As main result we obtained a closed 
formula for the Taylor coefficients for the diagrams under consideration, which allows their 
easy numerical evaluation with extremely high precision. The obtained coefficients were 
checked against those of Ref. 0, which were, however, only known explicitely up to a 
certain highest index in terms of rational numbers and £(2), £(3). 

This work was supported in part by the Volkswagen Stiftung under 1/71 293 and 
the BMBF under 05 7BI92P 9. We are grateful to M. Kalmykov, V. Smirnov and M. 
Tentyukov for discussions. 
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Figure 1. Planar diagrams with zero threshold occuring in Z —> bb and H - 
according to notation of Ref. || . Dashed lines are massless and solid lines massive. 
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